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Abstract 

Sum rules have played an important role in the development of many branches of physics since the 
earliest days of quantum mechanics. We present examples of one-dimensional quantum mechanical 
sum rules and apply them in two familiar systems, the infinite well and the single (^-function 
potential. These cases illustrate the different ways in which such sum rules can be realized, and 
the varying mathematical techniques by which they can be confirmed. Using the same methods, 
we also evaluate the second-order energy shifts arising from the introduction of a constant external 
field, namely the Stark effect. 

PACS numbers: 03.65.-w, 03.65.Db, 11.55.Hx 
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I. INTRODUCTION 



Quantum mechanical identities which relate time-dependent expectation values, influen- 
tial in the early days of quantum theory, continue to act as useful pedagogical tools in the 
modern curriculum. For example, the results often known as Ehrenfest's theorem(s), 1 

d(x) t d 2 (x) t /dV(x)\ 

can be used to show that time-dependent quantum expectation values are related to their 
corresponding classical equations of motion.- Identities restricted to time-independent ex- 
pectation values evaluated using energy eigenstates, \n), such as the quantum virial theorem 



(n\T\n) = In 
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and related hypervirial theorems,- are historically and pedagogically valuable as they too 
have clear classical analogs and can often be evaluated without resorting to direct integration. 

Similar relationships involving off-diagonal matrix elements, especially quantum mechan- 
ical sum rules, were also used to dramatic effect in the early days of quantum theory. For 
example, the Thomas-Reiche-Kuhn (TRK) energy-weighted sum rule,- again for energy 
eigenstates, 

Y^(E k -E n )\{n\x\k)\ 2 = — (TRK sum rule) (3) 
z — ' 2m 

k 

was used to describe the physics of electric-dipole interactions with atoms. It was originally 



obtained by requiring that the Kramers- Heisenb erg dispersion relation reduce to the Thomas 
scattering formula at high energies. Written in the form 

k k 

this was an important experimental check of the oscillator strengths (f n ,k) and an early 
confirmation of quantum results. Kramer was able to derive this relation in the context 
of matrix mechanics, reproducing the matrix version of the famous commutation relation 
[x, p] — ih.- Other early uses of sum rules included Bethe's study of energy loss mechanisms 
for charged particles in matter,- which made use of the relation 

fc2 2 

J2( E k ~ E n )\(n\e iqx \k)\ 2 = (Bethe sum rule) , (5) 
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eventually leading to the Bethe-Bloch formula. 

Since then, sum rules have been used in many areas of physics, including in atomic,- 
molecular, 8 solid state,- 1 ^ 1 ^ 1 ^ nuclear,— M.dkd!L and especially in particle physics.— iI&I2i22i2I 
One well-known paper applying sum rule methods to QCD22 is the 10 th most highly cited 
paper in the particle physics literature and over 2,000 papers on QCD sum rules have been 
published, with 60 appearing in 2007 alone.— 

The power of such sum rule identities is that they encode a large amount of information 
about both the energy spectrum and energy eigenfunctions of the system in a compact 
form, often in a way which is amenable to experimental confirmation. This in turn can 
probe assumptions about the fundamental interactions assumed or the calculation methods 
used to approximate physical systems. For example, QCD sum rules have been used to 
extract values of both the light and heavy quark masses, which are not otherwise directly 
measurable quantities. 

Despite their historical and contemporary importance, sum rules are not often treated 
in the context of standard quantum mechanics courses. The TRK sum rule is sometimes 
included in undergraduate quantum mechanics books,— but often only as a problem, and 
typically only using the harmonic oscillator. This lack of coverage may well be due to the 
paucity of tractable examples in familiar model systems to which students typically have 
exposure, or the level of mathematical analysis required to verify even the simplest cases. 

The purpose of this paper is to provide a suite of one- dimensional sum rules and to demon- 
strate the mathematical techniques required for their confirmation in two model quantum 
mechanical systems, the infinite well and the single (attractive) 5-function potential. In 
each case, the sum rules saturate differently and rely on different mathematical methods 
(summation techniques and contour integration methods) illustrating the diverse ways in 
which such sum rules are realized. The level of mathematical detail required, however, is 
low enough to be easily accessible to advanced undergraduate students. 

Explicitly confirming that these identities are indeed satisfied is not an empty exercise 
since it is possible to obtain surprising results, even from relatively simple systems such 
as the rigid rotator.— In addition, energy-weighted sum rule calculations are actually not 
exotic, since perturbation theory is discussed in standard textbooks in quantum mechanics. 
The expression for the second-order shift in the energy due to a perturbation V'(x) is given 
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by 

(2) ^ \{n\V\x)\k)\* 

k^n K 111 ^ ^k I 

which is a form of energy- weighted sum rule. Using this connection, we will find that we can 
make use of the exactly same techniques derived for confirming sum rules to evaluate the 
shift due to the addition of a constant external field, V'(x) = Fx, namely the Stark shift, 
in each of the model systems we consider. 

Introducing the concept of sum rules, using the tractable examples considered here, can 
certainly help students appreciate their use in later research applications. It can also help 
put the mathematical methods used in the same context as more familiar second-order 
perturbation theory calculations, one of the most important applications of undergraduate 
quantum mechanics, showing how the related required sums over intermediate states can be 
sometimes done in closed-form, and compared with explicit results. 



II. SUM RULE EXAMPLES 

The derivation of many energy weighted sum rules has been succinctly described^ as 
making use of a ' ...well-known technique which involves closure and evaluating a double 
commutator in two different ways.'' Such calculations rely on the fact that the solutions of 
the system under study form a complete set of states. For example, consider a system with 
energy eigenstates satisfying H\n) = E n \n). Then, for an arbitrary operator, O, we have 
the sum over off-diagonal matrix elements 

^|(n|6|A;)| 2 = ^(n|d|Jfe)(A;|d|7i) 

all k all k 

= (n\dlj2\ k )( k \\ 6 \ n ) 

{allk J 

= (n\d 2 \n). (7) 

We note that the sum over the complete set of intermediate states, \k), can include both an 
infinite sum (for discrete levels), an integral (for continuum states), or both. 

For the special case of O — x, we obtain the simplest dipole matrix element sum rule 
listed in Bethe and Jackiw,— >2I namely 

\J |(n|x|A;)| 2 = (n\x 2 \n) (x-closure sum rule) (8) 

k 
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with an identical sum rule relation for the off-diagonal matrix elements of the momentum 
operator. 

To derive the Thomas-Reiche-Kuhn (TRK) sum rule, we start with two commutation 
relations, 

[p, x] = — and [H, x] = [p 2 , x] = — p (9) 

% 2m mi 

where we assume a standard ID Hamiltonian of the form 

*=|Uk<«). (10) 
The first of these relations can be written in the form 

h 

— = (n\px — xp\n) = \ {(n\p\k) (k\x\n) — (n\x\k) (k\p\n)} (11) 

all k 

where we have inserted a complete set of states. The second relation in Eqn. ([9l) can be 



written as 



(n\p\k) = ^(n\[H, x]\k) = %m ^ En Ek \ n \x\k) (12) 
n n 



with a similar expression for (k\p\n). When used in Eqn. (JTTJ), this gives the desired result, 

h 2 



2m 



J2(Ek-E n )\(n\x\k)\ 2 . (13) 



k 

Wang 28 has derived a very general expression for the energy-difference weighted sum rules 
for the matrix elements of a well-behaved function of x, F(x), namely 
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Y J {E k -E n )\{n\F{x)\k)\ 2 = ^- { n 

k 



dF{x) dF\ 



x 



dx dx 



n ) (14) 



which simplifies if the function is Hermitian so that F(x) = F'(x). This general result can 
be used to immediately reproduce the TRK sum rule by using F(x) = x. We can then also 
derive the Bethe sum rule^ by using O = e iqx in which case we find 

J2(E k ~ E n )\(n\e^k}\ 2 = ^ ■ (15) 



k 

„2 



If we use F(x) = x , we obtain the so-called 'monopole sum rule,' which has been used in 
applications to nuclear collective excitations,— 

2h 2 

J2(E k - E n )\(n\x 2 \k)\ 2 = — (n\x 2 \n) . (16) 



k 



Wang 28 also discussed sum rules involving functions of the momentum operator, and 'mixed' 
x,p relations. 

Bethe and Jackiw 3 ^ 1 ^ derive several other sum rules for dipole moment matrix elements 
by using multiple commutation relations with the Hamiltonian, thus generalizing Eqn. (JSj), 
and yielding higher powers of the energy difference: 



y2(E k -E n ) 2 \( n \x\k)\ 2 = —(n\p 2 \n) 

k 

ih 2 

= —{E n -(n\V(x)\n)} . (17) 

(18) 



and 

Y^(E k -E n y\(n\x\k)\ 2 = — ln 

k \ 

where Eqns. f fT8|) and (fT9l) are described as the "force times momentum" and "force squared 1 
sum rules, respectively 

We note that not all of these sum rules are guaranteed to converge 3 ^ and in our case, 
because of the singular nature of the potentials used here as idealized models (the infinite 
well and the single-^ cases) several of these sum rules will not be applicable. 



E n f\{n\x 



fi 4 
2m 2 
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d 2 V{x) 
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dV(x) 
dx 



n 



(19) 



III. THE INFINITE SQUARE WELL 

The infinite square well (ISW) potential is the most frequently presented of all textbook 
examples of bound state systems and is frequently used as a model system to introduce 
students to tractable examples of research level physics, such as wave packet revivals.— We 
can confirm many of the sum rules discussed above for this case, making use of relatively 
straightforward mathematical techniques to evaluate the infinite sums which appear. (The 
only example we can find in the literature of the evaluation of sum rules in the context of 
the infinite well is a short discussion in Ref. [14].) 

We consider the standard ISW potential, defined by 

, for < x < a 
V(x) = { . (20) 

oo for x < and x > a 



The energy eigenstates and corresponding eigenvalues are 

, / x ■ fn-KX\ . _ h 2 n 2 n 2 . 

ip n (x) = \ - sin and E n = — — (21) 

V a \ a / 2mar 

where n = 1,2,... and the expectation value of x 2 required for the closure sum rule in 

Eqn. (jSD is easily calculated to be 

(n\x 2 \n) = a 2 (- . (22) 

X 1 1 1 \3 2n 2 n 2 J K ' 

The energy differences needed for the various sum rule calculations are given by 

while the off-diagonal matrix elements are given by 

2 f a . fmrx\ . fkirx\ 

\n\x\K) = — sin x sin ax 

a Jo \ a J V a J 



(24) 



k + n even 

-(8na/7r 2 )[A;/(A; 2 - n 2 ) 2 ] jfc + n odd 

so that for n even (odd) only odd (even) values of k will contribute. This result is due to 
the energy eigenfunctions' generalized parity property relative to the center of the well at 
x = a/2. For the closure identity in Eqn. (jSJ), we also need to include the diagonal matrix 
element, 

(n\x\n) = — . (25) 

This term does not contribute to the other sum rules, since the k = n term is suppressed 
by the (E^ — E n ) energy difference factor. In contrast to potential energy functions that 
are symmetric about the origin, such as the harmonic oscillator potential and the single 
5-function potential, the ISW potential as defined above is not symmetric and one must 
consider the k = n case for the closure identity. 

The position closure sum rule in Eqn. (JHJ) then reads 

allk v 7 k y ' 

where the summation is over even (odd) values k if n is odd (even). This is the first of many 
examples we will encounter where we require infinite summations of the form 

4 ±} w = E(p^ ( 27 ) 



where z takes on integral values, and where the summation is over odd, or even, 
values of k. For example, the required summation in Eqn. ( 1261) can be written in the form 



We provide a brief, but complete, review of how all of the sums required in this section can be 
evaluated using standard series expansions in Appendix |A] We note, however, that modern 
computer algebra systems (such as Mathematica) can easily handle such sums. Students 
may be allowed on first pass to use such tools and then asked to delve more deeply into the 
methods used to obtain the general mathematical results for this class of problems. 

For example, in modified Mathematica syntax, the summation over even integers k (rel- 
evant for n odd), yields: 

Sum [k~2/ (k~2-z~2) "4 , {k , 2 , Infinity , 2}] 

= (-12 Pi Cot [Pi z/2] - 6 Pi~2 z Csc[Pi z/2] "2 

+ 2 Pi~4 z"3 Cot [Pi z/2] "2 Csc[Pi z/2] "2 

+ Pi~4 z"3 Csc[Pi z/2] "4)/768z"5 

so that for odd integer values of z = n, we have (by hand or by using Assuming-> z G 
Integers in Mathematica) 



We obtain the same result (same function of n) for the summation over odd values of k 
(relevant for even n). A trivial modification (one character in fact) of the Mathematica code 
is all that is required. Using this result in Eqn. (I26I) . we then find that 




(28) 



k even 




(29) 



^\{n\x\k)\ 2 = \(n\x\n)\ 2 + J2\( n \ x \ k )\ 2 




(30) 



as expected. 

The TRK sum rule is then given by 




(31) 
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where the summation over k is only for even (odd) values for n odd (even). These sums can 
also be done in closed form and one finds 

k 2 



k even 
1 



(k 2 - Z 2 ) 



Si + \z) + z 2 St\z) 



\7T Z CSC , 

64z 3 L V 2 



TiZ 



2tt cot 



TiZ 



7l 3 Z 2 COt ( —J COS^ ^ — J 



TiZ 



(32) 
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A: odd 
1 



(Jfe 2 -3 



2A3 



7TZ 



7TZ\ 



n z z sec z ( — ) — 2ir tan ( — ) + tt 3 z 2 tan ( — ) sec^ ( — 



fKZ\ 



/TTZ 



64^ l" v 2 y - v 2 y 1 " ~ v 2 ; ~~ \2)\ (33) 

and we note the similarities in form. Inserting the appropriate odd and even values of n, in 
each case we find that 

T n + \n)=ll\n) = ^ 2 (34) 

for all integral values of n. This result, when substituted into Eqn. (I3T]) . directly confirms 
the TRK sum rule. 

Verification of the monopole sum rule in Eqn. (fT6l) requires a small, but important mod- 
ification of the summation methods. The off-diagonal matrix elements required for k ^ n 
are 

/ , 2,, x (~l) fc ~ n 8a 2 n / k \ 



7T" 



(k 2 — n 2 y 



while for k = n, one uses the result in Eqn. ( 1221) . Since the k = n term does not contribute to 
the sum (because of the associated energy difference factor) the left-hand side of Eqn. (|T6l) 
reduces to 

k 2 



E(^-^n)K 



n\x 



f h 2 ir 2 \ / 64n 2 a 
V2ma 2 J V vr 4 



E 



{k 2 — n 2 ) 3 



(36) 



k \ / \ / k^£ n 

and we must sum over all values of k ^ n since the even/odd pattern seen in the dipole 
matrix elements is not present in this case. 

In order to evaluate this sum, just as discussed in Appendix we can first generalize 
the sum to non-integer values of n, and then rewrite the sum as 



T(z:n) 



E 

k^n 



k 2 



(k 2 ~z 



2\3 



E 

all k 



k 2 



(k 2 -z 



2\3 



n 



(n 2 — z 2 ) 3 



(37) 



where the second term corresponds to the 'missing' term in the k ^ n summation. The first 
sum can be evaluated for arbitrary z, giving the result 

' TT COt(7T2) + 7T 2 Z CSC 2 (7TZ) — 2lT 3 Z 2 COt(7Tz) CSC 2 (tTz) \ 7l 2 



T(z;n) 



16z 3 



(n 2 — z 2 ) 3 



(38) 



Since we will be taking the limit where z — > n (and integral), we write z = n + e for general 
n, and we find that both terms have factors which diverge as 1/e 3 , 1/e 2 , and 1/e. If, however, 
we expand both terms about z — n (i.e. in small values of e) we find that these divergences 
cancel, leaving the finite result 

lim T(z; n) = lim T(n + e; n) = T(n) = Jj-j Q - (39) 

which when inserted into Eqn. (13"6T) reproduces the right-hand side of Eqn. (|16p . 

Many of the other sum rules discussed in Sec. [HI such as those that require derivatives 
of the potential energy function, Eqns. (fl8l) and (fl9l) . are not well-defined for the infinite 
square well (or the single 5-function in Sec. HVj) due to the singular nature of the potential 
energy function. While the matrix elements (n\e iqx \k) required for the Bethe sum rule in 
Eqn. ([SD are easily obtained in closed form, the summation methods discussed here are not 
immediately applicable. 

We can now use identical methods to evaluate the second-order shift of the energy levels 
of the infinite square well due to the addition of a linear potential, V'(x) = Fx, namely the 
Stark effect. In this geometry, where the ISW potential is not symmetric, the first-order 
energy shift is non- vanishing and is given by 

n F 1 

Eg* = (n\Fx\n) = . (40) 

The second-order shift has been evaluated by Mavromatis for the ground stated and then 
extended to a general state^ 1 ^ by using variations on the Dalgar no- Lewis method.— If we 
explicitly write the standard expression for the second-order energy shift, we have 

F ( 2 )_ y \(n\Fx\k)\ 2 / F 2 2ma 2 \ / 8na\' P 

which is formally identical to the class of summations discussed here. Using either Math- 
ematica or the results of Appendix |A] we find that the required sum (for either n even or 
odd) is given by 



Ek 2 _ 15n 2 n - % A n 3 



k 

so that 



(k 2 - n 2 f 3072n 7 

< (x / 15-(ri7r) 2 



1 H 2 { 247T 2 n 4 ' 
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The overall n-dependent form agrees with the results of Mavromatis,— £L22. who considered 
the related problem of the symmetric infinite well, for which the first-order correction van- 
ishes. This result is interesting in itself as the second order shift for the ground state is 
negative (as it always should be, since all states contributing to Eqn. (|41|) are higher in 
energy) but for n = 2 and higher, the shift changes sign. This is in contrast to the behavior 
of the harmonic oscillator, where the second-order shift is always negative, independent of 
quantum number. 

IV. THE SINGLE ^-FUNCTION POTENTIAL 

Another popular model system in which to investigate sum rule and perturbation theory 
results is the single (attractive) 5-function potential, defined here by 

V s {x) = -gS{x) . (44) 

The use of 5-function potentials as simply soluble models of potential barriers or wells has 
a long history in quantum mechanics, going back at least to Kronig and Penney 3 ^ who 
considered a ' series of equidistant rectangular barriers' and then took the limit where the 
' ...the breadth b of these barriers is made infinitely small and their height Vq infinitely large...' 
while not actually using the ^-function notation. 

Morse and Feshbach 35 explicitly considered the form in Eqn. (j44p . make note of the correct 
(dis) continuity condition on the energy eigenfunction at the origin, namely 

^(0+)-V/(0-) = -^(0), (45) 

cite it as being 'useful in the study of nuclear forces,' and go on to discuss the single bound 
state as well as scattering solutions. Frost 36 considered single and multiple attractive 5- 
function potentials as models of 1 hydrogen-like atoms,' the hydrogen molecule-ion and more 
complex systems. He was perhaps the first to explicitly comment on the similarities of the 
energy eigenvalue and eigenfunction for the single bound state of this system to the ground 
state of the Coulomb problem. Since then, single and multiple ^-function potentials have 
been widely used in model calculations in both the pedagogical and research literature. 37 

We note that compared to the two other most widely used simple ID models, the infinite 
well and harmonic oscillator, the 5-function potential has the advantage that it admits both 
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bound and continuum solutions, as does the Coulomb potential, and so it presents new- 
features compared to purely discrete spectra. 

The single bound (E < 0) state for the potential in Eqn. (I44p is given by 

-K \x\ 



i/> (x) = ^/K e~ Kolxl (46) 
where K = ^ with the corresponding bound state energy eigenvalue 

E °-~^--^- (47) 

One can then note many comparisons to the ground state of the hydrogen atom, if one 
defines the Coulomb potential as V c {r) = — ^ and one defines and substitutes a = in 
Eqns. (1461) and (14"T|) . Not only does the form of the ground state energy in Eqn. fj471) match 
that of the Coulomb potential, but the form of the energy eigenfunction in Eqn. (j4*6]) does 
as well. 

For use in confirming the closure relations in Eqn. ([8]), we find that for the ground state 
energy eigenfunction we have 

(0|x 2 |0) = ^. (48) 
The E > continuum states can be classified by their parity and are given by 

ipt\x) = -^sm(kx) (49) 



4 + \x) = 1 (K sm(k\x\)-kcos(kx)) (50) 

yir(k + K ) 



both of which have the same free-particle energy E^ = h 2 k 2 /2m. The combination of the 
single bound state in Eqn. (06]) and the continuum states in Eqns. ( T49l) and (1501) have been 
explicitly shown^ to form a complete set of states. The effect of the continuum states on a 
simple perturbation theory calculation has also been demonstrated by Kiang.— 

For the various sum rules, we will consider here the \n) = |0) case only, as others using 
purely continuum states do not converge. Because of the symmetry of the system, parity ar- 
guments dictate that the only non-zero dipole matrix elements connecting the single ground 
state to the continuum will arise from the ip^\x) states, and we find that 



(OMJfcH) = £J (V^e-^W) x (± Bin(te)) dx = 4^ JJff^ • ( 51 ) 
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The energy differences are then given by 

E k -E = ^-(k 2 + K 2 ). (52) 

and we note the similarities in form between these two expressions and the corresponding 
results for the ISW in Eqns. (|23|) and (|24l). 

The dipole matrix element closure relation in Eqn. ([8]) then becomes 

Ei<ow*<->>r=(^) f oct*^ ho|i2|0) (53) 

where the integral can be done by standard methods, and agrees with the value in Eqn. ( j48f ). 
The left hand-side of the TRK sum rule in Eqn. ([3]) gives 

r°° h 2 /ifi/r 3 \ r +OD k 2 h 2 

I iEt -E omx ^ )?dk =^(^)l 5^*4 (54) 

also as expected. We note the similarity in form of these integral expressions to the sum- 
mation results for the infinite well in Eqns. (l3Tj) and (T26|) . 

In order to confirm the monopole sum rule in Eqn. (Tl6]) . we require the off-diagonal matrix 
elements of x 2 for which only the even continuum states in Eqn. (|50|) contribute, giving 



:0|^I^)-8,/ t( /; a1) (55) 



We then find that 



E<*-*)i<oi»Wi'= g) P)f (^*=4 («) 

and we note that factors of (k 2 + Xq) from the energy difference and energy eigenfunction 
normalization in the numerator and denominator respectively cancel. 

The Stark effect for the single 5-function potential has been analyzed using exact results 
from the Airy function solutions,— 1 ^ as well as using the Dalgarno-Lewis method. 42 Using 
the dipole matrix elements derived above, we can evaluate the second-order energy shift 
directly, using the same kinds of straightforward integrals encountered so far. We find that 
l(2) f \(0\Fx\k^}\ 2 



fc( Er-4 o) ) 

2mF 2 flGK$\ f°° / 1 ] k 2 



dk 



& V / Jo I (k 2 + K 2 ) J {k 2 + K 2 Y 

(57) 
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which agrees with the results of Refs. [40]- [42], when put into this notation. We note that 
the entire contribution to the Stark shift in the ground state energy in this case comes from 
the continuum states and this result is one of the few examples of the explicit evaluation of 
the contribution of the continuum terms in such a calculation. 

We recall that for the hydrogen atom ground state, the total second-order shift 4 ^ can be 
written in the form 

^(H-atom) = -|^) (58) 

using the form of the Coulomb potential and perturbation theory, this result comes from 
summing over the contributions of both the bound states and continuum states. Ruffa 44 
has evaluated the continuum contribution to the expression in Eqn. ( 1581) in terms of a single 
integral and finds a net contribution of 0.4184 to the total 9/4 = 2.25 value of the pre- 
factor. It is perhaps then fairer to compare the second-order Stark result in Eqn. (I5l 
namely 5/8 = 0.625, to that partial contribution. 
Finally, the Bethe sum rule is given by 

k 2 q 2 



B = J2(E k -E )\(0\e^\k)\ 2 = -^- (59) 



2m 

k 

and in this case we will have two contributions to the left-hand-side, coming from the even 
(e) or odd (o) continuum states, namely 



B = B e + B 

(E k -E ) |(0|cos(gx)|fc (+) )| dk 



2 



+ / (E k -E )\(0\sm(qx)\k { - ) )\ dk (60) 
Jo 

and we consider each term separately. The first matrix element of interest is 

X = (0|mn(gaO|A; ( " ) ) 

e~ K ° x sin(gx) sin(kx) dx 
2kqK Q 



4JC, 



TV 







AK 



[(k + q)2 + K*}[(k-q) 2 + KZ}\ (61) 
where we use the symmetry of the energy eigenf unctions to evaluate the integral over positive 
values of x only. Recalling that E k — E = h 2 {k 2 + K^)j2m, we find 

(UK£\ r k 2 (k 2 + K 2 ) 

° 2m{nJJ [{k + q y + K 2 Y[{k- q y + K 2 y aK - {K)l} 
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Use of an integrated mathematics package (again Mathematica) returns the correct value for 
the integral, provided one correctly interprets the many cautionary restrictions on the values 
of K and q. Given the relatively complicated nature of the intermediate results coming 
from such programs, however, it is again important to be able to check the expressions 'by 
hand.' In this case, it simply involves extending the integral over the entire real line (since 
the integrand is an even function of k) and then using contour integration methods (see 
Appendix iBl for details), giving 

'KW\ \K 2 + q 2 /2\ 



b " ' 1 2m) \ K* + q 2 
For the even case, we require the matrix element 



(0|cos(gx)|£; {+) ) 



n(K 2 + k 2 ) 
tt(K 2 + k 2 ) 



e ° x cos(qx) [K sm(kx) — kcos(kx)] dx 



-2kK q 2 



[(k + q) 2 + K 2 ][(k - q) 2 + K 2 } 



(64) 



and the even contribution to the sum rule becomes 



{E k -E Q )m C os{qx)\k^)\ 2 dk 



K 2 q 2 \ (SKlq 2 



+oo 



A- 2 



dk. 



2m J V vr J [(k + q) 2 + K 2 } 2 [(k-q) 2 + K 2 } 2 

The integral can again be done with similar contour methods giving the result 

_ (h 2 q 2 \ f q 2 /2 
° e \2mJ\K 2 + q 2 

which can be combined with Eqn. (1631) to give 

H 2 q 2 



(65) 



(66) 



B = B + B e 



2m 



(67) 



as expected. 



V. CONCLUSIONS AND DISCUSSION 



We have presented an array of familiar (and not-so-familiar) one- dimensional sum rules, a 
number of which have proved useful in the development of many fields of physics. Using two 
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standard model systems as testbeds, we have illustrated the diverse ways in which such sum 
rules are confirmed, emphasizing the different mathematical techniques (infinite summation 
tricks and contour integration methods) used in each case. While the evaluation of the 
necessary summations or integrals can be simplified by the use of integrated mathematics 
programs, we have also provided the details necessary to demonstrate the same results from 
first principles. 

We have also noted the striking similarities of some of the expressions which arise for 
the same sum rules in the infinite square well and single ^-potential cases. Despite the 
qualitatively very different physical behavior of the two systems, they both begin with free- 
particle solutions. The infinite wall boundaries of the ISW force quantized eigenstates with 
E n = h 2 k^/2m, while the attractive 5-function gives the identical dispersion relation for the 
E > states, but with continuous fc-values. The 5-function case also includes one E < 
state for which the sinusoidal solution is analytically continued to the localized exponential 
form in Eqn. f j46l) . The connections between these two model systems are seldom if ever 
stressed, but appear very naturally in these sum rule calculations. 

We hope that the suite of exemplary problems discussed here can be useful to instructors 
in lectures as well as for homework problems, in both the advanced undergraduate and grad- 
uate quantum mechanics curriculum, especially by putting this important tool of theoretical 
physics into a historical and research context. 

APPENDIX A: INFINITE SUMS FOR THE SQUARE WELL PROBLEM 

Many of the sum rule and second-order perturbation theory results in Sec. II III for the 
infinite square well involve the evaluation of infinite sums of the forms 

s ^)=H w hy - ^M = £(^ < A1 > 

even k odd k ' 

where the both expressions are eventually evaluated using integral values of z = n, with n 
odd and even respectively so no divergences occur. While multi-purpose computer programs 
such as Mathematica can recognize and correctly evaluate such sums, it can be important 
for some students (and many instructors) to also be able to derive them 'from scratch.' To 
that end, in this Appendix we provide a very brief, but self-contained and complete, review 
of the mathematical tools necessary for their derivation from more basic results with which 



16 



students at this level should be quite familiar. 
We begin by considering the general expression 

OO .J 

s p& = E w^y P (A2) 

where the summation is over all positive integer values of k. The basic result we require is 
for the p = 1 case, namely 

" 2. (pTT) "2^ aT" ( ) 

fc=l v y 

which appears, for example, in Gradshteyn and Ryzhik.— This standard 'handbook' result 
can, in turn, be derived at a more fundamental level from a Fourier series expansion 46 by 
evaluating the Fourier components of the expansion 

OO 

cos(zx) = + [a n cos(nx) + b n sin(nx)] (A4) 

n=l 

over the interval (— tt, +vt); note that here z is considered a constant. The Fourier coefficients 
can be evaluated using standard integrals and we obtain 



, . simziT 
cosizx) = 



OO r 



E 



Z7T 

n=l 



2zsm(irz) cos(mr) 
7r(n 2 — z 2 ) 



cos(nx) (A5) 



since the b n = by symmetry. If we then specialize to x = it, and use the fact that 
cos 2 (mr) = 1, we find 

1 °° 1 

7 rcotM = --2^— — - (A6) 



n=l y 



and we note that this partial fraction expansion of cot(7r,2) correctly encodes the information 
on the divergences of the function at all integral (positive, negative, and zero) values of z. 
Rewriting this expression, we find that 



Si{z) " E (*nb) = h{l- n cot(7TZ) ) = i 

n=l v ' x ' 



7Tfot.(7r:) 



2z 2 2z 

confirming the handbook result of Eqn. (1A3j) . Such sums are already useful in that they can 
be used to evaluate quantities such as the Riemann zeta function, defined by 

oo oo ^ 2 

COO = E - g ivin § C(2) = Si(z = 0) = J2- = ^r (A8) 

n=l n=l 
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as the z — > limit of Eqn. (IA7I) . (Such results can be directly connected to integrals 
which appear frequently in the evaluation of quantities related to blackbody radiation which 
students encounter in standard textbooks 47 on statistical mechanics and can therefore be 
reinforced through such examples.) 

If we differentiate the result in Eqn. (1A7H with respect to z, we find that 

-S x {z) =2zJ2 J^T—Iy = 2zS ^ ( A9 ) 



so that in general 



1 dSJz 



s *«w = s^r (aio » 

thereby generating sums of arbitrarily high power. For example, this gives 

„ . . esc 2 (ttz) f-2 + 2ir 2 z 2 + 2 cos(2tt^) + ttz sin(27rz)l . . _ , 

S2(z) = '- > 11 (All) 

implying that C(4) = S 2 {z = 0) = tt 4 /90. 

Since our interest is often in summations restricted to the even or odd integers, we write 

OO ^ OO ^ OO 

*M - E (jiTpj - E + E (F^y - + s{-> W . (Ai2, 

k=l v 7 k even K ' k odd K 1 

We then note that 



(k 2 - z 2 ) ((2/) 2 - z 2 ) 

!y 1 Jo 

^(i 2 -(#) 4 1 

1 7TCOt(7T^/2) 



k even 
oo 



2^ 2 Az 



(A13) 



which then gives 
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^ / x -,f + i 7T r /vrz\ i 7r tan(7T2;/2) 
5} j (^) = S x {z) - Sy\z) = — cot — -2cot(7rz) = f — ^ (A14) 



42 



V 2 / 

where we use half- angle formulae in the last step. Both of the expressions in Eqns. (1A13|) 
and ( 1A14I) can, of course, be confirmed using Mathematica. 

The sums over higher powers of even/odd values of n required to evaluate Sp\x) and 
Sp~\x) in Eqn. (lAip are then obtained by repeated use of the differentiation trick in 
Eqn. (1A10I) . For example, we obtain results such as 

„(+),, CSC 2 (7T2/2)[-4 + tt 2 z 2 + 4cos(7T2;) + TTZSm(lTz)] 

S 2 (*) = ^ (A15) 
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and 



(_). . _ Tisec 2 (Txz/2)\Tiz -sin(7rz)] 



APPENDIX B: CONTOUR INTEGRALS 

The explicit evaluation of the integrals in Eqns. (|62|) and (1651) by contour integration 
techniques can be done by extending the region of integration over the entire real-axis. A 
contour consisting of a semi-circle of radius R can be then used as the integrands both have 
simple (double) poles at Zq = ±q + iK in the upper-half plane. The contribution to the 
contour integral over the circular arc vanishes as R — > oo, leaving 



/+oo 
F{k)dk = 2m^n i (Bl) 
-oo 



where the residues are given by 




^=tAtUt) [{z - zfYF{z)\) (B2) 



z (i) 



for = z^ and where in this case n 
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